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Abstract 

The general purpose spherical nonmagnetic detector ( SND ) is now taking data at 
VEPP-2M e~^e~ collider in BINP ( Novosibirsk ) in the centre of mass energy range 
of 0.2 -r- 1.4 GeV. The energy calibration of the Nal(Tl) calorimeter of the SND 
detector with cosmic muons is described. Using this method, the energy resolution 
of 5.5 % (fj) for 500 MeV photons was achieved. 



1 Introduction 



Electromagnetic scintillation calorimeters are an important part of many el- 
ementary particle detectors. One of the most complicated problems of such 
calorimeters is the energy calibration, i.e. determination of coefficients, needed 
for conversion of electrical signals from the calorimeter into corresponding en- 
ergy depositions, measured in units of MeV. Commonly used for that purpose 
are particles, producing known energy depositions in the calorimeter counters, 
e.g. gamma quanta from radioactive sources, cosmic muons, or particle beams. 

Cosmic radiation is a continuos and freely available source of charged particles. 
At the sea level the main part of them is represented by muons (~ 75%), 
having high enough energies (~ IGeV) to traverse the whole detector. In the 
calorimeter counters muons loose energy due to ionization of the medium. The 
average ionization losses are close to 5 MeV/cm in Nal(Tl). With reasonable 
accuracy one can also assume, that the trajectory of the muon is a straight 
line and that energy is left mainly in the counters crossed by this line. 
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Cosmic calibration of the SND calorimeter is a preliminary stage before the 
final calibration. Its goal is to obtain approximate values of the constants 
for the calculation of energy depositions. These constants are needed to level 
responses of all crystals in order to obtain equal first level trigger energy 
threshold over the whole calorimeter. They are also used as seed values for a 
precise OFF-LINE calibration procedure, based on analysis of Bhabha events. 
The requirements for the cosmic calibration procedure are the following: the 
calibration procedure must take not more than 12 hours, it must be indepen- 
dent of the other detector systems, and statistical errors in the conversion 
coefficients must be less than 1%. 

The scheme of this procedure was briefly described in [1]. Its important feature 
is that it makes use of virtually all muons detected in the calorimeter, greatly 
reducing the time required for calibration. In many other implementations of 
cosmic calibration, e.g. in the L3 detector at LEP [2], muons, traversing the 
detector only in certain directions, were selected. 



2 The SND calorimeter 

The SND detector [3,4] (Fig. 1) consists of two cyhndrical drift chambers for 
charge particle tracking, calorimeter, and a muon system. 

The main part of the detector is a three-layer spherical calorimeter, based on 
Nal(Tl) crystals. The pairs of counters of the two inner layers, with thickness of 
2.9 and 4.8 Xq respectively, where Xq = 2.6 cm is a radiation length, arc scaled 
in thin aluminum containers, fixed to an aluminum supporting hemisphere. 
Behind it, the third layer of Nal(Tl) crystals, 5.7 Xq thick, is placed (Fig. 2). 
The total thickness of the calorimeter for particles, originating from the center 
of the detector, is equal to 34.7 cm (13.4 Xq). The total number of counters 
is 1632, the number of crystals per layer varies from 520 to 560. The angular 
dimensions of the most of crystals are = = 9°, the total solid angle is 
90% of 47r. 

The electronics channel (Fig. 3) consists of: 

(1) phototriode with an average quantum efficiency of the photocathode for 
the Nal(Tl) emission spectrum of about 15% and gain of about 10 [5]. The 
light collection efficiency varies from 7 to 15 % for different calorimeter 
layers, 

(2) charge sensitive preamplifier (CSA) with a conversion coefficient of 0.7 V/pC, 

(3) 12-channel shaping amplifier (SHA) with a remote controlled gain that 
can be set to any value in the range from to a maximum with a reso- 
lution of 1/255. 
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(4) 24-channel 12 bit analog to digital converter (ADC) with a maximum 
input signal f/max = 2V, 

(5) in addition the SHA produces special signals for the first level trigger, the 
most important of which is the analog sum of all calorimeter channels — 
the total energy deposition. 

Each calorimeter channel can be tested using a precision computer controlled 
calibration generator. The amplitude of its signal can be set to any value 
from to 1 V with a resolution of 1/4096. The equivalent electronics noise of 
individual calorimeter counters lies within a range of 150-350 keV. 



3 Algorithm of the detector calibration with cosmic muons 



The calibration procedure is based on the comparison of experimental energy 
depositions in the calorimeter crystals for cosmic muons with Monte Carlo 
simulation. The simulation was carried out by means of the UNIM0D2 [6] 
code, developed for simulation of e^e~ colliding beam experiments. Simulated 
were muons crossing the sphere with a radius of 80 cm, centered at the beam 
crossing point, and containing all calorimeter crystals. The cosmic muon gen- 
eration program is based on data from [7]. Experimental data on the energy 
and angular spectra of cosmic muons were approximated using the following 
formula: 



dN 
dpde 



Fq ■ COS 9 ■ exp 



ln^(p/po) 
2cr2 



where Fq — 3 ■ 10~^cm~2s~''^ster~^(MeV/c)~"'^ — the flux of muons with a 
momentum oi p = Pq = 500McV/c at 9 = 0, 9 — angle between the muon 
momentum and vertical direction, k = 2.86, a = (2.86 — cos^^))/1.54. The 
data from [7], on the relative abundances of positive and negative muons, 
were approximated using the following expression: 



N+/N_ = 1,0 + 0,35 -exp 



-0, 22 In^ 



p 



7GeV/c 



(2) 



The passage of particles through the detector was simulated by means of 
the UNIM0D2 code. As a result the energy depositions in the calorimeter 
crystals were calculated in each event. A total of 1 000 000 cosmic events were 
simulated. 

The simulated events were processed as follows. The track of cosmic muon 
was fitted to the hit crystals and expected track lengths in each crystal were 
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calculated. Energy depositions and track lengths k were summed over event 
sample for each crystal separately and for each crystal their ratio was calcu- 
lated: Cmc(MeV/cm) = (E(MeV))/(l(cm)). Brackets denote averaging over 
event samples. Then the calorimeter channels were cahbrated electronically, 
using the generator. In order to account for the electronics nonlinearity, the 
calibration data were taken at two significantly different values of SHA gain 
and wide range of calibration generator amplitudes, covering the whole dy- 
namic range of the ADCs. In addition, data were taken at a working value of 
the SHA gain. 

Using all these data and assuming that the dependence of the ADC counts S'adc 
on the generator amplitude f/gcn is described by a second order polynomial, 
and that the generator amplitude itself linearly depends on generator code, 
one can obtain a, b, c, p, constants, defined by the following expressions: 

•S'adc = a[/gg„ + bUgen + C and Ugen = Ku(G + p), (3) 

where G and p arc a generator code and pedestal respectively, - the relative 
SHA gain. The next step is to obtain from these intermediate results the final 
constants A, B, and C for the working values of SHA gain, from the knowledge 
of which one can calculate the equivalent generator code G for any measured 
ADC count: 

G = ASl^, + BS^, + a (4) 

It relates an ADC count to a corresponding generator code, which allows the 
linearization of the ADC response. 

The next stage is processing of the experimental events. About 1.5 million 
events (Fig. 4) arc needed to be collected in special data taking runs with total 
duration of 4.5 hours. The processing itself is the same as that of the simulated 
events, but now the initial ADC counts are first transformed into equivalent 
generator codes. As a result, values of Cexp(Gen/cm) = (G(Gen))/(l(cm)) for 
each crystal are obtained. The meaning of C^xp is a code, which shoTild be 
written into the generator to produce signals on CSA input, equivalent to an 
average input signal produced by a 1 cm long muon track. 

Using Cexp (Gen/cm) and Cnic(MeV/cm), for each crystal we can calculate 

Ceai(MeV/Gen) = C^e(MeV/cm)/Ce.p(Gen/cm) - (5) 

an equivalent generator code, corresponding to an energy deposition of 1 MeV 
in a crystal. Then the gains of the SHA channels are adjusted to equalize 
contributions of different crystals into the total energy deposition signal and 
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the final generator calibration pass is carried out. It yields the coefficients 
Ok, 6k 5 and Ck, needed to transform S'adc into MeV according to the following 
formula: 



E{Mey) = akS^^, + bkSadc + Ck- 



(6) 



The use of the precision generator as a reference calibration source is very 
convenient not only because its linearity is much higher than that of the 
ADC, but it also simplifies the replacement of any element of the electronics 
channel, except the CPA or generator itself. All what is needed after that is 
to recalibrate the calorimeter with the generator. 



4 The description of the event processing algorithm 

The goal of the experimental and simulated events processing is to obtain the 
normalized ratios Cmc(MeV/cm) and Cexp (Gen/cm). The need for normaliza- 
tion of the energy deposition to the unit track length in a crystal arises from 
the fact that the raw energy deposition spectra in crystals are very wide and 
have very weakly pronounced maxima. On the other hand, such a normal- 
ization reduces systematic errors due to possible differences in the angular 
distributions of simulated and real muons, and differences in the directions 
of crystals axes over the calorimeter. The simulated and experimental events 
are processed using the same computer code. The only difference is that be- 
fore processing of a simulated event it is checked for passing the experimental 
trigger. The following selection criteria are used for event selection: the total 
number of hit crystals must be from 5 to 25. Crystals with energy depositions 
smaller than the threshold value (currently 5 MeV) are discarded. Similarly 
discarded are crystals with ADC counts less than 20. An event is selected if 
more than four hit crystals survive these cuts. For those crystals a straight 
line is fitted, using a least squares method, i. e. the sum of squares of distances 
from the crystal centers to the line is minimized. The line is parametrized as: 



where R is an arbitrary point on the line, v is a unitary vector in the direction 
of the line, t — parameter. Then, the squared distance from the center of i-th 
crystal to the line is: 



f^R + vt, 



(7) 



L^ = {X,-Ry-{v{X,-R)y, 
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where is a radius-vector of the crystal center. The minimized function is: 



F^^{R,v) = Y.LHal (9) 



where (7i is assumed to be equal to a half height of an i-th crystal. One point 
on a line can be determined immediately: 



^0 = TTTTZ^- (10) 



The direction of vector is determined by a maximum of the quadratic form 

F{Ro, v) = - ^o)) V^f ■ (11) 



In other words, v must have the same direction as an eigenvector of A, matrix 
of this quadratic form, corresponding to its maximal eigenvalue. Here 

^mn = E((^r - Ro)iX^ - Rl)/<^ll (12) 



m, n = 1, 2, 3 — coordinate indices. The equation for v was solved using an 
iterative procedure: 

v^+i^ Av^/\Av^\. (13) 



As seed value the unitary vector in the direction of the line between the centers 
of the two most distant crystals in the event is used. On each iteration step the 
presence of hit crystals at distances greater than 2a\ from the line is checked. 
If such crystals exist, two options are considered. The first one is to remove the 
most distant crystal from the list, the second one is to remove the crystal, most 

— * 

distant from the point Rq. For both cases the fitting of the line is repeated 
and the final choice is based on the minimal value of Finin(/2o, ^^)- If there are 
more than one crystal to be discarded then the event is rejected completely. 
Prom the remaining events only those are selected where F.^in{RQ,v) divided 
by the number of hit crystals is less than 0.7 (Fig. 6). 

After that, the expected lengths of the muon track in each crystal are 
calculated and averaged over the whole event sample together with the energy 
depositions E-^^. The values of average track lengths (ZJ and energy depositions 
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{Ei) are used then to calculate the energy deposition per unit track length 
Ci = {Ei) I (li) . Here i and j are the crystal and the event numbers, respectively. 

To estimate the statistical error in Q, the event sample is divided into groups 
of 50 events each and within each group the ratio Ri — (Ei) / (li) is calculated. 
Then, the error ctri in Ci can be estimated as a^i — {{Rf) — {Ri^y^'^/VN, 
where N is a number of groups. 

The total CPU time, required for processing of 1.5 • 10^ cosmic events is 2.5 
hours on VAXstation 4000/60. 



5 Events processing results. Compctrison of experimental and sim- 
ulated distributions 

The mean (Cmc) and their RMS values arc shown in the Table 1. The necessity 
of normalization of energy depositions in crystals to a unitary track length was 
studied. To this end the events with muons going through drift chamber were 
processed, and the corresponding C^^. coefficients were calculated. The ratios 
of (Cmc/ C^c) their RMS are also presented in the Table 1 together with 
the corresponding ratios of mean energy depositions {E/E^). The data show, 
that the normalization significantly reduces the dependence of coefficients on 
the angular distribution of muons. 

The distributions of experimental events over number of hit crystals with 
energy depositions higher than threshold value and over F^[^{Rq,v) are in 
good agreement with the simulated ones (Fig. 5, 6). Small differences between 
experimental and simulated distributions may be attributed to inaccuracies in 
the angular and energy spectra of the primary particles and first level trigger 
simulation. This also may cause small differences in the distributions over the 
angle relative to vertical direction (Fig. 7). Shown in Fig. 8 are distributions 
over total energy deposition in the detector, normalized to a unitary track 
length also agree well. 

A comparison of average track lengths of muons in the calorimeter crystals for 
the experimental and simulated events was carried out. The results are listed 
in Table 2. The statistical uncertainty of the Ccxp coefficients is less than 1 % 
for all three layers of the calorimeter. At the same time the uncertainties of the 
corresponding coefficients for simulated events are determined by simulation 
statistics. 
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6 Energy resolution of the calorimeter. Implementation of the cal- 
ibration procedure 



The energy resolution of the calorimeter was studied using e'^e~ — > e^e~ and 
e+e" — s> 77 processes. The energy distributions for electrons are depicted in 
Fig. 9. The events with polar angles of the particles in the range 45 < < 135 
degrees and acoUinearity angle less than 10 degrees were selected. To estimate 
the energy resolution quantitatively the spectra in Fig. 9 were fitted by a 
function 



.,.)....p|-i.[(!^)^,]} 



(14) 



?7 = 1 + 



t^{E - E^) _ sinh(tA/IK4) 
(7 t,t/\nA 



where is a particle energy, E'm — the position of maximum, A — normaliza- 
tion coefficient, tx - the asymmetry parameter, cr is a full width of the distri- 
bution at half maximum divided by 2.36. A, i^m, o", and tx were treated as free 
parameters during fitting. Energy resolution was defined as a-^/E = a/Ej^. 

Experimental and Monte Carlo resolutions of the calorimeter at 500 MeV are 
respectively 3.7 % and 5.5 % for photons and 3.5 % and 5.2 % for electrons. 
Peak and mean values of the experimental and simulated distributions agree 
to a level of about 1 %, but the widths of experimental distributions are signif- 
icantly larger. The possible explanation of the broadening of the experimental 
spectra could be electronics instability, systematic errors in calibration proce- 
dure and nonuniformity of light collection efficiency over the crystal volume. 

Study of the electronics and photodetector stability showed, that it can at- 
tribute to a maximum of 1 % shifts for the time duration of the collection of 
energy deposition spectra. Systematic biases were estimated by comparison 
of the calibration coefficients with those obtained using a special procedure 
based on minimization of the width of the energy deposition spectrum in the 
calorimeter for e'^e~ events. Such a calibration was carried out using experi- 
mental statistics collected in 1996 [8]. The average bias of coefficients relatively 
to those obtained with cosmic calibration was 3 % for the first layer. No sta- 
tistically significant bias was found in other two layers. The RMS difference 
in calibration coefficients, obtained using cosmic and e~^e~ calibration proce- 
dures is close to 3-^5 %. The resulting energy resolution after e+e" calibration 
was close to 5.0 % for photons and 4.6 % for electrons, which is still worse 
than that expected from Monte Carlo simulation. 

A Monte Carlo simulation was carried out taking into account a nonuniformity 
of light collection efficiency over the crystal volume. The results are shown in 
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Fig. 10 together with the experimental distribution. The energy resolution for 
simulated events decreased to 4.1 %. It shows that the broadening of experi- 
mental spectra was mainly determined by the nonuniformity of light collection 
efficiency over the crystal volume. 

The cosmic calibration procedure was used during a five month experiment 
with SND at VEPP-2M in 1996 [8]. Relative shifts in calibration coefficients 
are shown in Fig. 11. It can be seen, that for a one week period between 
consecutive calibrations, the mean shift of the coefficients is less than 2 % and 
RMS of their random spread is 1 2 %. 



7 Conclusion 

With the help of the described procedure of SND calorimeter calibration a 
statistical accuracy of 1 % in calibration coefficients and energy resolution 
of 5.5 % for 500 MeV electromagnetic showers was achieved. Mean energy 
depositions agree with simulation at a level better than 1 %. The total time, 
required for calorimeter calibration is not more than 5 hours. 
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Table 1 

Cmc coefficients. Comparison of calibration results based on different muon samples 
(MC). 



Layer number 


C'mc 


crc(MeV/cm) 


(C'mc/ C'ic) 


a 


{E/E') 


a 


I 


5.38 


0.14 


0.99 


0.02 


0.83 


0.04 


II 


5.43 


0.09 


0.96 


0.03 


0.71 


0.05 


III 


5.61 


0.08 


0.99 


0.03 


0.68 


0.05 



Table 2 

Comparison of mean track lengths of muons in calorimeter crystals for the experi- 
mental and simulated events. The statistical uncertainty of Cmc- 



Layer number 


((Zmc)/(/exp)) 


a 




I 


1.01 


0.02 


1.4% 


II 


1.00 


0.01 


1.0% 


III 


1.00 


0.01 


0.9% 
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Fig. 1. SND detector, section along the beams; 1 — beam pipe, 2 — drift chambers, 
3 — scintillation counters, 4 — light guides, 5 — PMTs, 6 — Nal(Tl) crystals, 7 — 
vacuum phototriodes, 8 — iron absorber, 9 — streamer tubes, 10 — 1cm iron plates, 
11 — scintillation counters, 12 and 13 — elements of collider magnetic system. 





Fig. 3. Electronics channel of the SND calorimeter: Nal — Nal(Tl) scintillator, 
VPT — vacuum phototriode, Preamp — charge sensitive preamplifier, Att — com- 
puter controlled attenuator, Shaper — shaping amplifier. Generator — calibration 
generator, ADC — analog to digital converter. Trigger — first level trigger. 
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Fig. 5. Distributions over number of hit crystals for the events with energy deposi- 
tions above certain threshold. N — number of hit crystals. Circles — experimental 
events, histogram — MC simulation. 
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- MC <F,„/N> = 0.38 
OEXP <Fm„/N> =0,39 
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Fig. 6. Distribution over residual values of minimized function in experimental and 
simulated events. 




Fig. 7. Distribution over the track angle with respect to a vertical direction. 
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Fig. 8. Total energy depositions in the calorimeter per unit track length in experi- 
mental and simulated events. 
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Fig. 9. Energy spectra for 500 MeV electrons. £^ is a measured energy, Eq is the 
beam energy. 
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Fig. 10. Energy spectra for 500 MeV electrons. The nonuniformity of light collection 
efficiency over the crystal volume was taken into account in Monte Carlo simulation. 
E is the electron measured energy, Eq is the beam energy. 




Fig. 11. The calibration coefficients spread. Points — average ratio of the current 
calibration results to the preceding ones, error bars — FWHM of the distributions 
of these ratios over the whole calorimeter layer divided by 2.36. Horizontal axis 
shows the number of days elapsed from the first calibration. Shown are the results 
for the second layer. Other calorimeter layers behave similarly. 
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